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,—1 The Hamilton action principle, also known as the principle of least action, and Lagrange equa- 
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CN tions are an integral part of intermediate and advanced undergraduate inechanics. Although the 

Q Hamilton principle is oft stated as "the action for any nearby trial world-line is greater than the 
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action for the classical world-line", the landscape of action in the space of world- lines is rarely 
explored. Here, for three common problems from introductory physics - a free particle, a uniformly 
accelerating particle, and a simple harmonic oscillator - we present families of trial world-lines, char- 
acterized by a few parameters, that evolve continuously from their respective classical world-lines. 
With explicit analytical expressions available for the action, they permit a graphical visualization 



ments the traditional Lagrange equation approach and is well-suited for advanced undergraduate 
students. 



I. INTRODUCTION 

Multivariate calculus and its applications to extremal problems are commonplace in un- 
dergraduate physics and mathematics.'^ A function /(xi, . . . , x„) of n variables is stationary 
at a point if all partial derivatives of that function vanish at that point: df/dxi = for 
i = 1, . . . ,n. Whether that point is a local maximum, minimum, or neither is determined 
by the symmetric second-derivative matrix Aij = d'^f/dxidxj = Aji at that point.'^ Thus, 
the problem of locating the stationary points of a function is equivalent to the problem of 
solving n equations in n variables, and can be tackled analytically or numerically. A com- 
plementary approach for locating the stationary points is to look at the global landscape 
of the function /(xi, . . . , x„) in M". This approach is typically emphasized^ only for n < 2 
when it is straightforward to obtain the graph (n = 1) or the surface plot {n = 2) of the 
function using ubiquitous software packages.^ Graphical exploration of the function in the 
vicinity of a stationary point also complements the standard technique used to determine 
whether the stationary point is a maximum, minimum, or a saddle point. However, this 
process becomes increasingly cumbersome when n > 3 and is therefore not emphasized in 
such cases. 

The Hamilton principle^ (also known as the principle of least actiorp'^ is an integral 
part of advanced undergraduate physics. This principle states that for a particle in one 
dimension, among all world-lines x(t) that connect two specified events {xi,ti) and {x2,t2), 
the Hamilton action 

S[x{t)] = f L{x{t),x{t),t)dt (1) 

is stationary for the classical world-line Xc{t) that obeys Newton's second law.'^^'^ Here, 
L{x,x,t) is the Lagrangian that, in general, depends upon the position x{t) and velocity 
X = dx/dt of the particle, and may also explicitly depend upon time. Since the action 5'[x(t)] 
is now a functional of the world-line x{t) - essentially a function of infinite variables x{t) for 
ti < t < t2 - the criterion for the stationary world-line Xc{t) is a differential equation, also 
known as Lagrange equation. 
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The Hamilton action principle for particles and fields in arbitrary dimensions provides an 
intuitive way to emphasize the deep connections^ between diverse topics such as mechan- 



ics, electromagnetism, and quantum mechanics where the requirement of a stationary action 
leads to Newton's second lawP Maxwell's equations^ and the time-dependent Schrodinger 
equatiorP respectively. There have been increasing efforts to introduce this topic in intro- 
ductory physics courses.'^'^ These efforts, based on reorganizing the introductory mechanics 
content, are supplemented by interactive graphical software that "hunts" for the classical 
world-line through sampling.'^ Recently, the question of whether the action for the classical 
world- line is a minimum or a saddle point, and the fact that it is never a maximum, have 
been addressed in detailP However, even in the simplest case of a particle in one dimension, 
complete numericaP^'^ and graphical exploration of the global landscape of the action in 
the space of all world-lines is impossible. Therefore, although the Hamilton principle is oft 
statecF'^ as "if you calculate the action for any world-line close to the classical world-line, 
it will be greater than the action for the classical world-line" , this exercise is seldom carried 
out explicitly. 

In this article, we present families of trial world-lines characterized by a few parameters 
for a particle in constant, linear, and quadratic potentials. The classical motion of the 
particle in these cases is familiar to advanced physics undergraduate students. In each case, 
the trial world-lines in a family are continuously connected to the corresponding classical 
world-line. Using these trial world-lines, we graphically and analytically explore the action 
landscape in the vicinity of the classical world-line. In addition to explicitly verifying that 
the action is minimum for the classical world-line, this exercise permits a look at corrections 
to the action that arise when the trial world-lines deviate from the classical one. Due to the 
global nature of this approach, we are able to explore situations in which there are multiple 
classical world-lines that connect the two given events. We emphasize that by construction, 
the families of trial world-lines presented here are only a subset of the set of all nearby trial 
world-lines, and thus the action landscape obtained here is not exhaustive. This exercise 
provides, however, a way to explore the action landscape graphically using common software 
programs such as Mathematica, Grapher (Mac), or Sage.l^ 

II. THE ACTION LANDSCAPE: FROM A FUNCTIONAL TO A FUNCTION 

In this section we consider three common problems from introductory physics: a (almost) 
free particle, a uniformly accelerating particle, and a simple harmonic oscillator .'HUSl In each 



of these cases, the classical trajectory of the particle, typically obtained by using Newton's 
laws, is well-known to the students starting an intermediate or advanced mechanics course; 
therefore, in the following, we do not present the derivation of the classical world-lines. 

A. Free particle 

Let us start with the case of a free particle of mass m in one dimension. Since this system is 
translationally invariant, without loss of generality we choose (xi,ti) = (0,0). The classical 
world-line connecting the initial event (0,0) and the final event (^2,^2) is given by Xc{t) = 
X2t/t2, and, since the Lagrangian for this system is L{x,x) = mx^/2, the corresponding 
free-particle action is given by Spp = mx\/2t2. Let us consider a family of trial world-lines 
Xa{t) = X2{t/t2)°' characterized by a single parameter a. It is straightforward to evaluate 
the action as a function of a 



S{a) = Spp 



a^ 



(3) 
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Note that a well-defined action requires a > 1/2. The action is minimum when dS/da = 
SFp2a{a — l)/{2a — 1)^ = or, equivalently, a = 1. Fig. [I] shows typical trial world-lines 
along with the classical world-line at a = 1 (solid line), and dimensionless action S{a)/SFP 
as a function of the exponent a. We see that the graph has a single minimum at a = 1. 
A Taylor expansion of the action (|3| shows that for nearby world-lines, a ~ 1, the change 
in the action 6S = S{a) — Spp varies quadratically with the distance from the classical 
world-line, 6S = Spp^a — 1)^, although, as we will see below, this conclusion is not generic. 

Another family of trial world-lines for the same problem is given by x^(t) = X2 sin{ujt)/ sin{ujt2) 
for cj 7^ and Xuj=oit) = X2t/t2. Note that Xuj(t) is an even, continuous function of u. The 
Hamilton action for the trial world-line is given by 



S{uj) = S_ 
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The trial action diverges at ut2 = nn, and has a single minimum at nn < ujt2 < {n + l)7r 
between consecutive divergences. It follows from Eq.Q that the global minimum occurs 
when ujt2 G (— 7r,7r). Since the action is an even function of the parameter u, it follows 
that it is either a minimum or a maximum, but not a saddle point, at w = 0, the classical 
world- line. Taylor expansion of Eq.Q around a; = shows that the change in the action 



SS for this family of trial world-lines is quartic in the distance from the classical world- 
line, 5S = +Spp{ut2y/4:5. Figure [2] shows typical trial world-lines, along with the classical 
world- line a; = (solid line), and the dimensionless action S{u)/Sfp as a function of the 
dimensionless parameter ut2. We see that, indeed, the action S{uj) has a very flat minimum 
at w = 0. These two families of trial world-lines show that even for a free particle, the action 
landscape in the vicinity of the classical world-line is very rich. 

B. Free particle with a perfectly reflecting barrier 

Now, we consider the free particle with a world-line that starts at {xi,ti) = (0,0) and 
ends to the right at X2 > at time ^2 > 0, in the presence of a perfectly reflecting barrier at 
Xh < 0. In this case, the particle is free at all times except when it collides with the barrier; 
at that instance, its momentum is reversed while its kinetic energy is unchanged. This 
system provides the simplest example where there are two classical world-lines that satisfy 
the given boundary conditions, and without further information (say, about the direction of 
the initial velocity) it is impossible to choose a unique classical world-line. It is useful, for 
the purposes of analytical calculations, to define the barrier position as Xb = —^7x2 where 
r] > is a dimensionless scale factor. Due to the presence of the barrier, in addition to the 
classical world-line Xci{t) = X2t/t2, there is another world-line that represents the particle 
initially moving to the left, bouncing off the barrier, and then traveling to the right, 

,,, / -VX2t/iCt2) 0<t<et2, ,^, 

Xc2{t) = < (5) 

[ -r]X2 + (1 + 7])x2{t - it2)/ [(1 - 0^2] it2<t< t2. 

Here ^ = i]/{2ri + 1) < 1/2 represents the fraction of the time that the particle spends 
moving to the left. Note that the speed of the particle along the second classical world- 
line Xc2{t) is [21] + 1) times larger than its speed along the first classical world-line Xci{t). 
Hence the actions for the two classical world-lines are given by Si = mx\/2t2 = Spp and 
5*2 = (2?7 + 1)^5*1 respectively. Since we have already studied the action landscape around 
the first world-line, we will focus on the second. To explore the action landscape in the 
vicinity of the second classical world-line, we consider a family of trial world-lines that are 
parametrized by two exponents 

^PnV^) = \ (6) 

-r]X2 + (1 + r/)a;2 [{t - ^^2)7(1 - 0^2]'' ^^2 < t < ^2- 



Note that, as before, a well-defined action requires /3,7 > 1/2, and that the second classical 
world-line corresponds to (/3,7) = (1,1)- The Hamilton's action for the trial world-line is 
given by 



S{P,l) = S2 



P' . ,. .^ !■ 
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^7T77p^ + (l-0 



(7) 



(2/3-1) ^ "^(27-1) 
This action is a sum of actions for the two paths characterized by /3 and 7 weighted by the 
fraction of the time that the particle spends on each path. Since the trial action is now a 
function of two variables, it is stationary when dS/dfi = = dS/d^. It is straightforward 
to verify that (/3,7) = (1, 1), which corresponds to the second classical world-line Xc2(t), is 
the only solution. Since the time that the particle spends moving left cannot be greater 
than the time it spends moving right, ^ < 1/2 < (1 — ^), the action S'(/3,7) is anisotropic 
around (/3,7) = (1, 1), and the anisotropy vanishes only when ^ — t- 1/2. Figure [s] shows a 
typical contour plot of the dimensionless action S{f3,'y)/S2 in the (/3,7) plane for a particle 
with barrier at X^ = —8x2/4 or, equivalently, C, = 0.3. As expected, we see that the action 
is minimum at (/3,7) = (1, 1). 

Next, we can consider a set of trial world-lines that span both classical world-lines: 
Xp,a,i3,'y{t) = pXa{t) + (1 — p)xfs^^{t). Notc that the first classical world-line Xci{t) is ob- 
tained when {p,a) = (1, 1) whereas the second classical world-line Xc2(t) is obtained when 
(p, (3, 7) = (0, 1, 1). We leave it as an exercise for the Reader to obtain the action S{p, a, /3, 7): 
in order to simplify the calculation, choose a = /3 = 7; plot the dimensionless action 
S{p,a)/ Spp in the (p, a) plane and verify that the action has local minima at the two 
classical world-lines. 

C. Uniformly accelerating particle 

We now consider a particle of mass m in one-dimension being acted upon by a constant 
force F. Since this system, too, is translationally invariant, we choose (xi,ti) = (0,0). The 
classical world-line joining the two events is given by 

where A = Ft^/2mx2 is the dimensionless acceleration and u = t/t2 G [0, 1] is the dimension- 
less time variable. Since the Lagrangian for this system is L{x, x) = mx'^/2 + Fx, the action 
for the classical world-hne is given by Sp = {mxl/2t2) [1 + 2yl - A^/?,] = Spp [1 + 2A- A'^/3] . 



It increases continuously from its free-particle value when < A < 3, and is monotoni- 
cally decreasing for A > 3 becoming negative when A > A^. = (3 + 2a/3). Motivated by 
Eq. (Is]) , we consider a family of trial world- lines characterized by a single parameter 7 and 
the corresponding action, 

x^{t) = X2 [(1 - A)u^ + Au^^] , (9) 



^(7) = s. 
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+ T-. ^ + ^ 7^ + V7^ 7T + 



(10) 



(27 -1) (47 -1) (37 -1) (27 + 1) (7 + 1) 

Note that once again, a well-defined action requires 7 > 1/2. The action is stationay for 
any value of A when dS/d'-f = or equivalently when 7 = 1, the classical world-line. A 



Taylor expansion of Eq. ( 10 ) near 7 = 1 shows that for nearby trial world-lines, the change 
in the action is positive, SS = +Sfp{20A'^/27 - A + 1)(7 - 1)^ > for any A, and thus 
confirms that 7 = 1 is a minimum. Figure |4] shows the trial world-lines for A = 1 and A = 7 
including the classical world-lines (solid lines), and the dimensionless action S{'~f)/\SF\ as a 
function of 7 for different values of A. The action has a single minimum at 7 = 1, and the 
minimum value changes from +1 to —1 when A > Ac as is expected. Thus, Eq.([9]) provides 
a suitable one-parameter family of trial world-lines to explore the action landscape for a 
uniformly accelerated particle. 

For the same system, an alternate family of trial world-lines parameterized by two expo- 
nents and the corresponding action are given by 

XaAt) = ^2 [(1 - ^)m" + AU^] , (11) 



S{a,(i) = S, 



FP 



[1-Afa^ ^2/32 2A{1-A)al3 AA"^ AA{1 - A) 



(2a -1) (2/3-1) (a + /3-l) (/3 + 1) (a + 1) 



.(12) 



Note that, as expected, Eq.(|12j) reduces to Eq.(lO) when (3 = 2a and that a well-defined 
action requires a,(3 > 1/2. The trial world-line Xa^pit) is continuously connected to the 
classical world-line, Eq.([8]), that occurs when (a,/?) = (1,2). The action is stationary when 
dS/da = = dS/df3. It is possible, albeit tedious, to analytically verify that (a, (3) = (1, 2) 
is a solution of these equations, and that the matrix of second derivatives at this point has 
positive eigenvalues confirming that it is a minimum. Figure [5] shows a typical contour plot 
of the dimensionless action S{a, P)/\Sf\ in the a-/3 plane for A = 3 < Ac- The action 
shows a clear minimum at the classical world-line (a, /3) = (1,2). We also note that action 
landscape in the vicinity of the classical world-line becomes increasingly anisotropic as A 



increases. Thus Eqs.([9]) and (11 ) provide two distinct ways to graphically explore the action 



landscape in the vicinity of the classical world-line, Eq.Q. 

D. Simple harmonic oscillator 

As the last example, let us consider the particle of mass m in a simple harmonic potential 
given by V{x) = mojQx'^/2. The classical world-line Xc(t) connecting events (xi,ti) and 
(x2,t2) is given h^^^^ 

..w=x. "''"°'r''' +^. """°^''' (13) 

sm uJo[t2 - tij sm Uo{t2 - ti) 



Note that Eq.(13) is invalid when iX'o(t2 — ti) = rnr or, equivalently when X2 = ±xi. In these 
two cases the action for the classical world-line is zero and thus it is impossible to quantify "a 
small change to the classical action" without resorting to the fundamental (quantum) unit 
of action h. Since the Lagrangian for this system is given by L{x,x) = mx'^/2 — mu}Qx'^/2, 
it is straightforward to evaluate the action for the classical world-line and we gel!^^'^ 

SsHO = 'muo [{xl + a;2)cosu;o(t2 - ti) - 2x1X2] /2sinu;o(t2 - ti), (14) 

= ,„(.; + x^H^A«[5|^], (15) 

where At = (t2 — ti), and Q = uoAt ^ rnr and B = 2x1X2/ {xl + X2) are dimensionless 
quantities. Since |-B| < 1 the classical action changes from positive to negative when Q > 



ilc = arccos(i?). Note that Eq.(13) satisfies the classical world-line boundary conditions for 



any u. Therefore we consider a one-parameter family of trial world-lines 

sma;(t2 — ti) sma;(i:2 — ti) 



sin2r2'u sinfiw 

2Qu Qu 



The corresponding Hamilton action is given by 

S{u) = ^(^i+^2KA^ 1(^2 _ ^) (^ -BcosQu) + {u' + 1) 
4 sin \lu { 

^ ' (17) 

where u = uj/uq is the dimensionless parameter. It is straightforward to analytically check 
that the action has a minimum at m = 1, dS/du\u=i = and (PS/du'^\u=i > 0. Figure p^ 
shows the typical dimensionless action S{u)/\Ssho\ for Q = 1.45 when B = 0.1 (solid line) 
and B = 0.2 (dotted line). In each case, the action has the minimum at m = 1. We recall 
that as B increases, Qc = arccos(i?) decreases; therefore, for a fixed value of Q, the action 

8 



minimuiii changes from +1 at B = 0.1 (solid line) to —1 at B = 0.2 (dotted line). We see 
from the graph that the action for a trial world-line deviates quadratically with the distance 
from the classical world-line, SS oc {u — 1)^; in addition, its curvature at m = 1 diverges as 
Q ^ Qc from both above and below. It is straightforward to extend this analysis to a family 
of trial world-lines characterized by two parameters, 

sin/i(t-ti) smu{t2-t) 

Xf,,u[t - X2- T\+Xl- 77 -T-. [i-O 

smfi{t2~ti) smz/(t2 — ti) 

In this family, the classical world-line occurs when (/i, u) = {uq, uq). We leave it as an exercise 
for the Reader to obtain an analytical expression for the action S{fi, v) and to graphically 
explore the dimensionless action S{ii,v)/\Ssho\ in the (/x, z/) plane in the vicinity of the 
classical world-line. 



III. CONCLUSIONS 

The principle of least action and Lagrange equations are two integral elements of interme- 
diate and advanced undergraduate mechanics. Subsequently, they are taught in advanced 
electromagnetism and quantum mechanics courses to emphasize the deep connections among 
the aforementioned subfields of physics. In all instances, however, an explicit calculation of 
the Hamilton's action for "nearby trial paths" is seldom carried out. In this article, we have 
graphically and analytically explored the action landscape in the vicinity of the classical 
world-line. Since the action is a functional of a given world-line x{t) that obeys the bound- 
ary conditions, it is impossible to exhaustively span all trial world-lines close to the true 
classical world-line Xc{t). Instead, for three familiar problems from introductory physics, we 
have presented families of trial world-lines that are characterized by a few parameters and 
that evolve continuously from their respective classical world-lines. Within each of these 
families, the Hamilton's action is reduced from a functional to a function of the (few) pa- 
rameters that characterize the family. The calculation of action for these trial world-lines 
is accessible to advanced undergraduate students. Its graphical exploration using common 
software^nSl permits an explicit, albeit limited, demonstration of the principle of least action 
with an approach that is complementary to the Lagrange equations.'^nSl 
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FIG. 1. (color online) (a) The left panel shows trial world-lines for a free particle, x^it) = x^itjt'i)^ ^ 
for a=0.6 (dotted line), q=1.8 (dashed line), and the classical world-line q=1 (solid line). Note 
that for trial world-lines with a < 1 the velocity diverges as t — )• 0, for trial world-lines with 
a > 1 the velocity vanishes as t — )• 0, and the action for the trial world-line is finite only when 
a > 1/2. (b) The right panel shows that the dimensionless action S{a)/SFP has the minimum at 
the classical world-line a=l. As expected, the minimum action is equal to that for a free particle, 
Spp = mx\/2t2, and has a quadratic expansion around the minimum, 5S = Sppia — 1)^. 
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FIG. 2. (color online) (a) The left panel shows trial world-lines for a free particle, x^{t) = 
X2 sin wi/ sin a;i2 for ujt2=1.6 (dashed line), ^^2=0.8 (dotted line), and the classical world-line 
ujt2=0 (solid line). Note that x^{t) is an even function of lo. Since the global minimum of Eq.Q 
occurs in the range — tt < u}t2 < vr, we only focus on |a;t2| < vr. (b) The right panel shows the 
dimensionless action S{u;)/Sfp restricted to a; > since S{uj) is an even function of w. The action 
shows the minimum at the classical world-line w=0, and the flatness of the curve near a; = is 
consistent with the quartic expansion, 6S{uj) = S'irp (a;t2)^/45. 
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Plot of S(y6,y) with clasical world-line at (y6,y)=(l,l) 
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FIG. 3. (color online) Typical contour plot of the action S{P,'j)/S2, Eq.Q, with trial world-lines 
a;^,y(t) for a free particle bouncing off a barrier located at X^ = —3x2/4 or, equivalently, ^ = 0.3. 
The classical world-line Xc2{t), Eq.([5|, corresponds to (/?, 7) = (1, 1) where the action shows a clear 
minimum. Note that as ^ — )• 1/2, the action anisotropy in the vicinity of the classical world-line in 
the /3-7 plane vanishes. 
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FIG. 4. (color online) (a) The left panel shows trial world-lines for a uniformly accelerated particle 
x^{i) = X2 [(1 — A)u^ + ^u^'''] where 7 = 1 corresponds to the classical world-line for any A = 
Ft\/2mx2- Note that as A increases the time the particle spends in the region x < increases; 
hence the action for the trial path decreases and becomes negative when A> Ac = {?> + 2v3). (b) 
The right panel shows the dimensionless action as a function of the trial world-line parameter 7 for 
various values of A. For A < Ac (dotted line and solid line), we see that the action has a minimum 
value of -|-1 at 7 = 1, whereas for A > Ac (dashed line), the action has a minimum value of —1 at 
7 = 1. In each case, the action increases quadratically as 7 deviates from the classical world- line 
value 7 = 1. 



14 



3.0 



2.5 



2.0 



1.5 



1.0 



Plof S(a,yg) with classical world-line at (a,/3)={l,2) 




FIG. 5. (color online) Typical contour plot of the action S{a, I3)/\Sfp\, Eq.(12), with trial world 



lines Xajsit) for a uniformly accelerated particle with yl = 3 < Ac- The action shows a clear 
minimum at the classical world-line (a,/3) = (1,2). We note that as the value of A is increased, 
the anisotrophy in the Q-/3 plane in the vicinity of the classical world- line increases. 
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FIG. 6. (color online) Typical plots of the action S{uj) for a simple harmonic oscillator for trial 



world- lines X(^(t), Eq(16). The action is plotted in units of absolute value of the classical action 
SsHO and the frequency uj is in units of the natural frequency ujq of the system. For all values of 
B = 2x1X2/(2;^ + X2) and ft = woAt, the action ^(cj) shows the minimum at a; = ujq. For a fixed 
0, = 1.45 as B is increased from B = 0.1 (solid line) to B = 0.2 (dotted line) , Qc = arccos(i?) is 
reduced, and hence the action minimum changes its value from +1 (solid line) to —1 (dotted line). 
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